extern doublereal epslon_(doublereal *x);
